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(g) Proc©d6 et dispositrf permettant a un modem de so synchroniser sur un transmetteur de donnees 
nume piques par voie hertzienne en presence de brouifleurs. 



@ Le proced* consiste a estimer (9) la matrice 
de correlation ftxx des signaux re$us du trans- 
rnetteur sur un ensemble de n capteurs de 
reception, a calculer (10) la matrice de correla- 
tion inverse ft& f a calculer (11) des vecteurs 
d'intercorrfeiation fxdi entre les signaux recus 
sur ('ensemble des n capteurs et un signal de 
replique d(k) connu. a calculer un critere de 
synchronisation (15) en effectuant les produits 
des vecteurs d'tntercorrelation transposes 
eonjugues ftii et de la matrice des vecteurs 
dlntercorrelauon ?xdi. et A comparer (16) la 
valeur du critere obtenu a une valeur de seuii i\ 
determines pour placer la synchronisation sur 
I'echantlHon pour lequel la valeur des criteres 
depasse la valeur de seufl ^. Pour estimer le 
decalage Doppler, la Sequence de synchronisa- 
tion est* decoupee en tranches d'un nombre 
determine de symboles pour effectuer des cai- 
culs de transformation de Fourier rapide sur 
des vecteurs d'intercorreLation determines sur 
chaque tranche. 
Application : Transmission radio de donnees. 
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La presents invention cencerne un procede et un disposrtif permectant a un modem de se synchroniser 
sur un transmetteur de donnees numeriques par voie hertzienne notamment en presence de brouilleurs. 

Dans les systemes de transmissions de donn&es numeriques par voie hertzienne la synchronisation est 
effective, a partlr de sequences de symboles connus du recepteur et intercalees entre les symboles d'infor- 
5 mation, par correlation entre le signal reeu et un signal de reference construtt a partir des symboles connus. Le 
resultat est compare d un seuil pour didder de la synchronisation. Des systemes permettant la synchronisation 
d'un modem a partir d'un seul capteur et d'un seul recepteur sont par exemple connus de Karticle D.V, Sarwate et 
M.B. Pursefy intitule "Cross correlation properties of pseudo random and related sequences" paru dans IEEE Vol. 
66, pp. 593-619, May 1980. Mais les techniques mises en oeuvre qui permettent en outre une estimation des ca- 
10 ract6ristiques du canal de transmission deviennent rapidement inefficaces en presence de brouillage. 

Lorsque le rapport signal a bruit n'est pas suff isant pour la prise de synchronisation par une technique 
mono-capteur. des techniques visant a accroitre le rapport signal d bruit et a permettre la synchronisation sont 
utilisees. 

' Ces techniques peuvent etre regroup&es en deux families. Une premifere famille vise a pretraiter le signal 
is recu avant d'effectuer la synchronisation du signal pretraite. La seconde famille vise a ameliorer le rapport 
signal d bruit et d effectuer la synchronisation conjointement 

Dans les techniques relevant de la premiere famille il faut ggalement distinguer deux categories : 
une premiere cat6gorie met en oeuvre un pretraitemeru uhJi^^ptoUr^ de brouilleurs 

avant d'effectuer la synchronisation. Ce traitement sst connu par exemple de I'anicJe de L.B. MILSTEIM 
20 "Interference rejection Techniques in spread spectrum communications", paru dans Proa IEEE, Vol. 76, 

pp. 657-661, June 1988. Suivant cet article on cherche a predire la contribution du brouilleur, dont la 
bande est supposee plus petite que la bande du signal, sans pr6dtre celle du signal utile. Cela permet 
d'amGllorer les performances de la synchronisation en presence d*interf6rences bande 6troite (type 
FSK). Mais le traitement reste inefficace en presence «fun brouillage occupant toute la bande utilis6e 
25 par le signal utile. Cela est en partte due au fait que les techniques de pnMraitement mono-capteur ne 

permettent un traitement des signaux que dans un espace limits a Taxe des temps et a Taxe des fre- 
quences. 

- La seconde categorie concerne les techniques mettant en oeuvre un pnMraitement multi-capteurs de 
rejection, precedent la prise de synchronisation. L'id6e consiste a mettre en oeuvre un f iltrage spatial, 

30 discriminant le signal et les brouilleurs par leur direction d*arriv6e, le temps, le spectre ou la puissance 

comme cela est decrit dans Kouvrage de RA. Monzingo, T.W. Miller, intituld "Introduction to adaptive 
arrays", edit© par John Wiley and Sons, New- York, 1980 et la these de P. Chevalier, intitul6e "Antenna 
adaptative : d'une structure liniaire a une structure non lineaire de Vol terra", Th&se de doctors t Uni- 
verstte de Paris Sud, Juin 1991 puis de correler le signal de reference avec la sortie de Tantenne. Les 
35 techniques de f iltrage spatial discriminant le signal et les brouilleurs par le temps, le spectre ou la puis- 

sance sont general ement sous-opti males, alors que celles exploitant la connaissance de la direction 
d'arrivee du signal necessitent la mise en oeuvre d'une technique de goniam6trie, coOteuse et delicate 
dans la gamme des frequences HF notamment du fait des multi-trajets. 
Pour ces rafeons, les techniques de la seconde famille, visant a effect uer la rejection de brouilleurs conjoin* 
tement avec la prise de synchronisation, a partir de la connaissance du signal de reference, apparaissent plus 
prometteuses. Une telle technique a par exemple 6te propos£e par M. Compton dans son article "An adaptive 
array in a spread spectrum communication system'' paru dans Proc. IEEE, vol. 66. pp. 289-298, March 1978, 
pour la synchronisation de signaux PN. a partir d'une architecture analog ique integrant des "Delay Lock Loop". 
Une seconde technique f ut proposee par L.E. Brennan et t.S. Reed dans leur article intitule "An adaptive 
45 array signal processing algorithm for communications", paru dans IEEE Trans- Aero-SysL vol. A£S~18, n°1. 
pp. 12^-130, June 1982 pour les transmissions numeriques etcompletee parD.M. Dugloset R.A. Schottzdans 
I'artide intitule "Acquisition of spread spectrum signals by an adaptive array", paru dans IEEE f rans-Acous- 
Speech-Signat-Proo, vol. ASSP-37, n"8, pp. 1253-1270, August 1989. Ces techniques apparaissent toutefois 
beaucoup plus cooteuses que les techniques mono-capteur et sont malheureusement inapplicables telles quet- 
&> les en presence des derives Doppler du signal inherentes aux transmissions ionospheriques. 
Le but de ('invention est de pa liter les in convenient s precites. 

A cet effet, ('invention a pour objet un proc£d6 permettant k un modem de se synchroniser sur un trans- 
metteur de signaux numeriques par voie hertzienne en presence de brouilleurs consistant a: 

- estimerla matrice de correlation Rxx de signaux re$us sur un ensemble de N capteurs de reception par 
55 blocs composes d'un nombre d*6chantillons d'une sequence de synchronisation transmise par le trans- 
metteur, ^ 

- a calcuier la matrice de correlation inverse Ri<\ 

2 
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- a calculer des vecteurs d'intercorrelation entre les signaux X(k) reeus surl'ensernble des N capieurs 
et un signal da replique d(k) connu. 

- a calculer un critere de synchronisation effectuant le produit scalaire entre la vecteur d'intercorrelation 
Txa et ie vecteur obtenu par produit de IMnverse de la matrice de correlation R« et du vecteur d'inter- 

5 . correlation r^, 

- a comparer la valeur du critere obtenu a une valeur de seuil r| determine© pour placer la synchronisation 
sur I'echantiilon pour lequel la valeur des criteres d 6 passe la valeur de seuil tj, caracterise en ce qu'il 
consiste pour estimer le decalage Doppler inherent aux communications entre emetteurs et r6cepteurs, 
a d^couper la sequence de synchronisation en un nombre determine de M tranches de ko symboles, a 

10 effectuer sur ces tranches des calculs de vecteurs d'intercorrelation partielles f XdJ (i = 1, M) de telle 

sorte que sur chacune des tranches de ko symboles, le decalage Doppler fasse tourner ta phase d'une 
valeur infferieure a 90°. a effectuer des calculs de transforrnee de Fourier rapide (FFT) sur les vecteurs 
fxai obtenus, et a calculer une valeur du critere de synchronisation a partir des vecteurs r m {kAf) t obtenus 
en resultat des calculs de transforrnee de Fourier rapide, et de la matrice de correlation Rx* 
15 (.'invention a egalement pour objet un dispositif pour la mise en oeuvre du procede. 

Le procede et le dispositif selon rinvention ont pour avantages qu'ils permettent a la fois d'assurer la syn- 
■" chrtjnisation cTun modem dans des conditions de brouillage severe et en presence de mulrj-trajets, d'estimer 
Ha^6ftvcen frequence ertsrer ernetteur et rScepteur et de deteetsrtes instants^*YMe-d^ 
■ * associes au signal utile, pourvu que-ceux-ci scient de puissance suff isamment forte par rapport a la puissance 
2o d e bruit de fond. 

D'autres caracteristiques et avantages de rinvention apparaTtront dans la description qui suit faite au re- 
gard des figures des dessins annexes qui represented : 

- La figure 1 un schema fonctionnel d'une antenne adaptative selon rinvention. 

- La figure 2 un exempfe de mise en oeuvre d'un filtre adapts spatial a discrimination par la modulation 
25 pour la realisation d'une antenne adaptative selon rinvention. 

* La figure 3 un schema d'antenne adaptative pour !a mise en oeuvre du precede de synchronisation spa- 
tiale selon {'invention. 

- La figure 4 un graphe pour illustrer la methode d'estimation de la matrice de correlation de signaux pro- 
venant de N capteurs. 

30 - La figure 5 un graphe pour illustrer le decoupage des symboles de ta sequence de synchronisation par 

tranche pour le calcul de vecteur d'intercorrelation en presence de decalage Doppler. 

- La figure 6 un graphe montrant revolution de critere de synchronisation mis en oeuvre par I* invent ion 
en presence <Tun decalage Doppler. 

- La figure 7 un organigramme retracant Tensemble des etapes nScessaires a la mise en oeuvre du pro- 
35 c6de de synchronisation selon rinvention. 

Le procede de synchronisation mis en oeuvre par ('invention fait appel aux techniques connues de syn- 
chronisation mono-capteur pour les appliquer a la realisation d'antennes adaptatrves sur replique. 

Selon la technique de synchronisation mono-capteur, le critere de synchronisation C(n) est obtenu en cal- 
culant la correlation y(n) entre K symboles connus dans une sequence de synchronisation et les K echantillons 
40 qui suivent I'arrivee d'un 6chantillon n sur le capteur, en normalisant les resul tats par la puissance p(n) du signal 
re$u. Ces calculs sont effectues en appliquant les relations : 

i K 

« y(n) = ~£x(k+n)cr(k) (1) 

1 K 

P(n)-^2x(k-Hn)x*(k + n) (2) 

C ( n)=i^ « 

55 ou: 

- x(k) est le signal issu du capteur utilise. 

- d(k) est le signal de replique connu (K symboles places en debut de trame). 
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- 9 represente V operation de conjugaison sur (es nombres complexes. 

Ces calculs sont reiteres pour chaque echantillon n, etlorsque le entire C(n) depasse un certain seuil ti, 
dependant de la probability de fausse alarme, la synchronisation est decidee. L'instant de synchronisation est 
determine pari© rang d'echantillon n 0 pour lequei C(Ho)> x\. La normalisation par la puissance p{n) permet de 
5 choisir un seuil independent du niveau de signal repj. 

Dans ce precede la synchronisation est obtenue lorsque le rapport signal a bruit est suff isanL Ses per- 
formances dependent en partie du nombre d'echantillons pris en compte pour la correlation. Cependant en 
presence de brouilleurs forts, la synchronisation ne peut etre effectuee. 

Uidee mise en oeuvre dans ('invention consiste alors a ef fectuer la synchronisation, en sortie d'un traite- 
io ment d'antenne adaptative. permettant la rejection des brouilleurs. 

Le schema fonctionnel d'une antenna adaptative est represente a fa figure 1 , celle-ci comporte un reseau 
1 deNcapteurs^ ... 1 N , un reseau def iltres adaptes 2 (constitues par des circuits de ponderations amplitude- 
phase, des f iltres a reponse impulsionnellefinieRIFou infinie RJ I) couple h un processeurdecommande adap- 
tatif 3. Le process eur 3 calcule les coefficients W n a W N du reseau de f litres 2 pour Padapter au signal utile 
1$ et aux brouilleurs, en I'absence de toute connaissance a prion sur les brouilleurs at a partir d'une information 
mrnimale sur le signal utile. Les coefficients W 1 a W N permettant de former un trou dans le diagramme de 
rayonnement du reseau d'antenne dans la direction des brouilleurs, tout en maintenant un certain gain dans 

-""-"Tdniirection^du signal utile.- ' - • — - . - — — = --«~* . — - - - 

Differents types d'inf ormation sur le signal utile sont exploiter par I'antenne adaptative. Ceux-d sont d'or- 
20 dre spatial lorsque la direction d'arrivee du signal utile est connue, d'ordre temporelle lorsque les instants de 
presence et d'absence de signal utile sont connus par exemple en evasion de frequence, reiatifs a la forme 
d'onde lorsqu'une information sur la modulation du signal utile est disponible pour disposer d'une repllque du 
signal utile en intercalant par exemple des sequences connues dans la trame utile, cf ordre spectral lorsque 
le signal utile occupe une certaine bande de frequence et que les brouilleurs occupent une bande plus impor- 
25 tante. ou statistiques pour separer les dtf ferentes sources arrivant sur le reseau en exploitanl par exemple les 
statistiques cTordre 4 des signaux en supposant que les differents signaux sont statistiquement independants. 

Dans le cas de I'antenne adaptative representee a la figure 2, rinformatton sur le signal utile est constitute 
par les sequences de sym boles connus intercalees dans la trame entre les sequences de symbol es dlnf or- 
mation. Les f iltres places en sortie de chacun des capteurs 1 1 & 1 n sont formes respecttvement par des circuits 
30 de ponderation amplitude-phase Wi a W N des signaux x,(t) a x^t) fournis par les capteurs. Les signaux obtenus 
en sortie des circuits de ponderation sont appliques a un circuit de summation 4. 

Les caracteristiques du reseau de f iltres 2 sont calculees en minirntsant une Erreur Quadratique Moyenne 
(£QM) estimee £, entre un signal de replique d(k) et le signal y(k) sortant de I'antenne de la fecon schematisee 
a la figure 2 ou les elements hornologues & ceux de la figure 1 sont reperes avec les memes references, % est 
35 def inie par la relation 

i = Zlv 00 - <*f = j Rx(k) - 4kf (4) 

« k=l k=f 1 

ou: "*"*■-. 

- y(k) est la sortie de I'antenne adaptative et est def inie par 

y(k) = WX(k) (5) 

- W est le vecteur poids definissant I'antenne adaptative, W* etant le vecteur transpose conjugue corres- 
pondanL (le signe t en exposant def init I'operateur correspondant) 

- X(k)'est le vecteur fbrm6 par les signaux re$us sur les capteurs, 

- d(k) est le signal de replique connu. 

La difference y(k) - d(k) est obtenue sur la figure 2 a la sortie d'un circuit soustracteur 5. 
Le vecteur poids W est calcule par la formule connue de Wiener : 

W = R^fxd (6) 

avec: 

Rxx "4 £x(k).xt(k) (7) 
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matrice de correlation des signaux provenant des capteurs 1 1 & 1 N . (X^k) etant le vecteur transpose conjugue 
du vecteur X(k)). 

1 K 

F Xd = ^£X(k)d(k) 
K k=l 

vecteur d'intercorreiation entre X(t) et d(t) (8). 

L'antenne adaptative ainsi utilisee conduit dune maximisation du rapport signal/(bruit + interference) 
ccnnu sous la designation SNIR abrdviation anglo-saxonne de "Signal to Noise + Interference Ratio 0 obtenu 
en sortie de l'antenne, avec une possibiltte de rejection pouvant aller jusqu'a N-1 brouilleurs, et est appelee 
Filtre Adapt© Spatial (FAS). 

La recherche de la synchronisation est effectuee selon 1'invention en sortie de l'antenne adaptative en 
discriminant le signal selon sa modulation, et en uttlisant le critere de synchronisation mono-capteur. Le sche- 
ma de principe de la synchronisation spatiale est represents 3 la figure 3 ou les elements homologues a ceux 
des figures 1 et 2 sont representes avec les memes references. Le entire de synchronisation devient en uti- 
*lSantfea ho^^ " " . - ^ 

C(n)4Sy(^n)d(k) (9) 

25 y(k**-nji represente la sortie de l'antenne adaptative conduisant £ une minimisation de I'erreur quadratique 

moyenne £(n) estim6e sur K echantillons entre la rtpfique d(k) et le signal arrivant pris a partir de t'6chantillon 
n ( suivant la relation 

K 

30 5(n) = 2|y(k + n) - d(k)| 2 avec y(k + n) - W(n)*. X(k + n) (1 0) 

k=1 

avec: A 
35 W(n) = Rxx- 1 (")-^(n) (11) 

Les differentes correlations sont calcul6es sur les K echantillons qui suivant rarrivee d'un echantillon n 
suivant les relations : 



to 



15 



40 



45 



55 



1 K 

Rxx(n) = ^r Sx(k + n)X l (k+n) (12) 
F Xd(n)=^i;x(k + n)d7k) (13) 



Lorsque I'echantiilon n correspond k la position de synchronisation, les poids W{n) def inissent le filtre 
so adapte, et y(k+n) contient alors le signal utile debarrasse de la contribution des brouilleurs. Dans le cas contrai- 
re, aucun signal n'est correle avec la replique d(k) et y(k+n) ne contient plus alors que la contribution du bruit 
de fond. 

^expression du critere de synchronisation spatiale se simplifie en ecrivant : 



1 K ♦ . 4 r. k ^ 

C (") = -j7 I>(n) l X(k+n)d (k)= W(n)* 1 £X(k + n) d"(k) 
*k«1 l K k-1 ) 



(14) 
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soit 

C(n) = ^(nJRxxtnJ-^n) (15) 

De meme que pourle cas mono-capteur. ce calcul est reiter$ pour cha que echantillon n. et lorsque le critere 
C(n) depasse un certain seuil q, dependant de la probability de fausse aJarme, la synchronisation est decides 
(la normalisation par la puissance est id obtenue par rintermediaire de Ri 1 )- 

L'instant de synchronisation est alors determine par le rang de l'6chanttllon pour lequel C(n 0 ) est superieur 
A une valeur de seuil determine* r\ K 

Les calculs correspondant sont effectues sur la figure 3 par le bloc 7 de calcul du critere C(n) et par le 
dStecteur de seuil 8, £ partir de la sequence de synchronisation d(t) et du signal y(t) correles par le correlateur 
6. 

Cela necessite pour chaque echantillon. d'une part de calculer N correlations entre le signal de repligue 
d(k) et les signaux x(k) regus par les N capteurs (soit N fois les calculs effectues pour la synchronisation mono- 
capteur). et d'autre part d*estimer puis d'inverserla matrice dMntercorrelation des signaux capteurs. 

Au cours de ces calculs la synchronisation et le calcul du filtre adapts 2 sont effectues conjointement 
Lorsque la valeur du critere depasse le seuil ti, calcule, la synchronisation est decidee et le vecteur poids de- 
f inissant le filtre 2 sen deduit facilement par la formule connue de Wiener. Le vecteur poids W peut ensuite 
^Vullil^pauf-filfcrfcr^ - 

En supposant que la replique d(t).est de puissance 1, les correlations tendent vers les valeurs suivantes, 
lorsque le n ombre <f echantillons K devient tres grand ; 

Rxx = R + n» S S 1 

Ou: 

- R est la matrice de correlation des interferences et du bruit de fond, 

- S est la matrice du vecteur dtrecteur du signal utile, 

- n, est Ea puissance du signal utile recu sur un capteur. 

Le critere de synchronisation qui converge alors vers une valeur Co s'ecrit encore: 

c ^StR-i S _ SNIRopt (m 

0 " 1 + SJR-i S ~ 1 + SNIRopt K ' 
ou SNIRopt est le SNIR asymptotique obtenu en sortie du filtre spatial adapt© 2. 

Lorsque le SM1R est fort, le critere tend vers 1 T lorsqu'il est faible, il tend vers 0. L'importance du choix du 
seuil ri, fixe la probability de detection de la synchronisation. Le seuil depend natureUement de la probability 
de fausse alarms souhaitee- a 

Le calcul du critere de synchronisation C. de r^ et de Rxx est effect u6 pour chaque position de la syn- 
chronisation et notamment la matrice Rxx est estimee et inversee egalement pour chaque position de la syn- 
chronisation, la puissance de calcul demandee est etev6e_ Mats il est possible de red u ire la puissance de calcul 
en estimant la matrice de correlation Rxx moins souvent sur un intervane de temps plus important (K' Echan- 
tillons au lieu de K avec IC>K de la fagon representee d ia figure 4) et en utilisant cette estimation de Rxx sur 
toutes les positions de synchronisation incluses dans cat intervalle de temps. 

Dans ces conditions la matrice Rxx s'ecrit : 



+ e (17) 



ou : 

- B(tfc) est la contribution du bruit et des brouilleurs sur X(t fc ), 

- S{t^ est la contribution du signal utile sur X(tO ? 

- Ko est ('instant d'apparition du signal utile sur le palter, 

- c est un terme prenarit en compte I'intercorrelation entre le signal B(t) et le signal utile S(t). qui tend vers 
0 lorsque K devient grand. 

Pour la position de synchronisation testae ; 



K* K 1 

XB(t k )B t (t k )+ SSOk)^) 



k=1 



k=K B 
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rxd — 



K 0 +K 

Z X(tk)d"00 
k=Ko 



(18) 



Dans les conditions asymptotiques (avec K»1), lorsque le signal utile est synchronise avec la re pi i que 



10 



15 



20 



25 



30 



d(t), Rxx et r Xd s© simplif ient et s'Scrivent : 



A 



R + + (19) 



f„ = -v/^S (20) 

Le critere obtenu en estimant la matrice de correlation Rxx sur K' echantillons possede comme valeur 
asymptotique, apres un calcul utilisant le lemme dlnversion matricielle: 

c . o = SNIRofit (21) 

j| ^ iC-Ko + 1 sNIFtopt 

K' — Kft f 1 

C omm r ^ V? , _ . < 1 % ajprs .CV j> Qc ' ^ 

La valeur du entire apparait dans ces conditions augmentee d'un facteur dependant de I'instant d'appa- 
rition du signal utile sur la frequence consideree, ce qui entrame une amelioration de la prise de synchronisa- 
tion, au prix d'une puissance de calcul reduite. 

D'autre part, le nombre K* ne doit pas etre choisi trap grand, sa valeur est determines pour que sur 
K* echantillons, renvironnement de brouillage puisse etre consider* comme stationnaire. 

Le seuil d prendre en consideration depend de la probability de fausse alarme (Pfa). La (Pfa) est la pro- 
babilite qu'un bruit soit detecte comme une synchronisation. Le calcul de la valeur du seuil n a lieu en calculant 
Pfa=P(C>n.) lorsque Tantenne recoit uniquement du bruit 

En presence de bruit le critere de synchronisation s'ecrit 



C = 



^ Z B k B k* I ^ZBk^ (22) 



35 



En supposant que B* est un bruit gaussien complexe de puissance a 2 , independant entre les capteurs, 
lorsque K est suffisamment grand, le terme du milieu de la relation (22) devient : 



40 



r 1 k "\ 

l K k=l 



= a +2 Id 



(23) 



45 



(Id etant la matrice identite ou unitaire) 

te critere de synchronisation peut alors etre approxime de la facon suivante : 




K 2 - 2 



£Bk<*c 
KM- J 



(24) 



so 



Soit: 



ss 



N 



C = 



Z B k ' ^ 



k=1 



N 



K 2 a 2 



Z(^+lf) US) 
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avec: 



5 r^^-^+jii- (26) 

&l represent© le bruit obtenu sur le capteur i, pour r^chantillon k. 

+ jl, peut etre approximg par un bruit gaussien complexe de variance K a 2 . Dans ces conditions R, et l| 
to sont chacun des bruits gaussiens rtels de variance 



75 



(K» 2 ) 



En not?.n»T{ ~ R h ou !; r C : 



20 



C "ST?|M (27) 



oii T t est une variable aleatoire qui suit une loi normale centime de variance La variable aleatoire Z x • 



C devient 



Ti suit done une loi nonmale centree de variance 1. 



30 



2N 



CS 2RZ2f (28) 



i=1 



35 



Comme la suite des 2| peut etre consideree comme une suite de 2N variables aleatoires independantes 
gaussiennes de moyennes nulles et de variances i , la variable aleatoire 



i=i 



4S 



suit une loi du x 2 * 2N degres de liberie. 
La probability d'avoir C S n s'ecrit alors ; 



p(C ^ n) = p 



so 



(2N - 

£zf *2Kt| 
M=1 



= Pfa (29) 



Le tableau suivant donne differentes valeurs de seuils r\ obtenues par recherche dans une table dormant 
P(Z=iA) pour une loi du x 2 £ 2N degres de liberty et pour Pfa = 1 CH. 



55 
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N = S 




N = 4 




N = 3 


s 


N = 2 



=> A = 42 => r| = 0,26 

^> A = 37 => n"0.23 

A = 33 ^ T| = 0,21 

=> A = 28 => ti^O.18 

Corome entre remission et la reception, le signal transmis peut Stre affect6 d'un certain decalage en fre- 
quence, ou decalage DoppJer, du esse ntielle men t aux derives en frequences des oscillateurs locaux des 
10 synthetiseurs a remission et d la reception, ainsi qu'aux conditions de propagation ionospherique, il est ne- 
cessaire de compenser ce decalage. Le risque est s'il n'est pas compense, d'empecher toute prise de syn- 
chronisation. 

En appelant Af 0 , te decalage Doppler entre le signal re$u par t'antenne et le signal attendu et en eonsiderant 
la position de synchronisation, le vecteur signal re$u peut s'ecrire avec les notations des paragraphes sui- 
75 vants : 

X(t) = ^ d(t) e*^' S + B(t) 
La rn a trice de g^f^atioq estiin6e pe ut aJo rs s'ecrire : ^ ^ ^ ^ ' - 

20 K ■ 

RXX S[V^d(t k )e 2 ^- Af oA S + BlfK^^ddk^^Vs + Bdk)] 1 (30) 



*XX = * + 1 £[tts|d(t k )| 2 S + B(t k ) B{t k ) t ] (31) 



ou £ est une matrice prenant en compte les intercorrelations entre ia contribution du signal utile et le bruit b(t). 
Lorsque K»1, dans les conditions asyrnptotiques, R^x peut s'approximer par : 

Rjoc =r «■ S S* + a 2 Id (32) 

La relation (32) montre que le decalage Doppler ne modifie pas le catcul de Rxx- En effet. le decalage 
s'elimine naturellement lorsque le signal utile est multiple parsa quantity conjuguee. 

Le vecteur d'intercorrelation peut s'ecrire : 

K 



? Xd = £ d(y^-*Mc S +B(tk)] cf (33) 

Fxd = «' + £ I^^JdCtkjj 2 or*J**M« s] (34) 



oil t' est un temne prenant en compte I'intercorrelation entre le bruit B(t) et la replique. Lorsque K»1 , rx 8 peut' 
s'approximer par : 



S (35) 



La relation (34) montre que le decalage Doppler modifie Testimation du vecteur d'intercorrelation r M et 
que du fait que la phase du terme en exponentiel complexe peut tourner de plusieurs fois 360° sur rintervaiie 
de temps utilise pour le calcul de la correlation, integration ne peut plus avoir Heu de fapon conerente ce qui 
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empeche le critere de synchronisation multicapteurs de converger. Le processus suivant permet d'estimer le 
decalage Doppler et de le compenser pour rendre possible la prise de synchronisation. 

Seion ce processus, les K sym boles de la sequence de synchronisation sont decoupes de la facon repre- 
sentee a la figure 5 en M tranches de Ko symboles (M Ko = K), de telle sorte que sur chacune des tranches 
de Ko symboles. le decalage Doppler fasse taurnerla phase du signal re$u par Tantenne <fune valeur inferieure 
a 90*. af in de permettre une integration coherente de r^ sur chacune des tranches. 

Si r X di designe la valeur de Tintercorrelation estimee sur la tranche i, celle-ci s'ecrit : 



r Xdi = 



^ k=iK 0 vl 



S i = 0 (36) 



En mettant en facteur dans la relation (36) un terme dependant uniquement de Ko, la longueur de la tranche 
elementaire consideree, et du decalage Doppler, la^re liation (36) devient : 



r Xdi = 



^i:>/^d((k + iKo)T)| 2 e^ F e 



e r - S i = 0,...,M-i (37) 



Ce qui petit encore s'ecrire ; 



r X di^x,ea*^ ! JK,.si=0,...,M-1 (38) 



F. 



La relation (38) montre que les vecteurs sont des vecteurs, echantillortnes a Leur transformee de 
Fourier discrete est alors : 



M-1 



s^**** FC/K ° (39) 
i=0 



2 j£^=£.IK 0 
>Xd(Af)s £ xje Fe ^ S 

i=0 



(40) 



La relation (40) montre que lorsque Af = Af 0 . le decalage Doppler est compense. ce qui rend maximum la 
valeur du critere de synchronisation r Xd (Af) t Rxx" 1 rxa(A0. 

Af in d'estimer Af, c'est-a-dire af in cf estimer le decalage Doppler qui conduit a la valeur maximal© du cri- 
tere, et afin de reduire le coQt de calcul correspondant, une transformee de Fourier rapide est effectuee sur 
les vecteurs (le nombre de tranches eJementaires M etant choisi pour former une puissance de 2) ; ce calcul 
donne: ' 



M-1 



At, 



-.ik M-1 



M k = 0 M-1 (41) 



oti af, est la resolution Doppler telle que : 

Af, = 



Fe/K 0 
M 



(42) 



Si la resolution est insuf f isante, une amelioration d'un facteur 2 peut elre obtenue en calculant une trans- 
formee de Fourier rapide. sur les vecteurs r Ml "decalees" en frequence de la facon suivante : 
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( ■ ilO 



-2 - 

e "M k = 0 M-1 (43) 



L'operation consists & multiplier chacun des vecteurs d^ntercorrelation r^ parle terme e-0 ,ri/M > et a effec- 
tuer ensuite un calcul de transformee rapide sur les vecteurs ainsi obtenus. 

Le critere de synchronisation est ensuite estime surchacune des cases de calcul des transformees de 
10 Fourier precedences par la relation : 

C ( k f)=^f)^^f) lc.OL_.2M-1 (44) 

Bten evtdemment, si la resolution est toujours insuffisante, il est possible de I'ameliorer en utilisant le 
* meme processus pour calcul er des vectaurs- ■ - — ■ t 



20 



oO n est un entier quelconque. 

25 Un test sur la case de la transformee correspond ant au maximum du critere de synchronisation est effec- 

tive. Si le seuil est depasse pour una valeur ko ^~ don nee du decalage Doppler, il est alors possible d'af finer 
le calcul du decalage Doppler (indispensable pour la demodulation) en calcuiant la valeur du critere pour dif- 

30 ferentes valeurs de decalage Doppler entourant Ko En effet, a Tissue de cette etape, la resolution sur Tes- 

timation du decalage Doppler est Pour rameliorer, il faut procederpar dictiotomie : a chaque iteration et 

de la facon representee a la figure 6 t la valeur du critere pour la valeur de decalage Doppler situee au milieu 
des deux valeurs de decalage Doppler dormant les valeura maxima du critere lors de Iteration precedente 
35 est calculee. La recherche par dtchotcmie est arretee des que la resolution souhaitee sur 1'estimation du de- 
calage Doppler est obtenue. Apres n iterations de Talgorithme de recherche par dichotomie, la resolution, qui 

etait initialement de ~ % est : 

«ill = Jk- ( 45) 

L'enchaTnement des stapes de calcul precedente pour la mise en oeuvre du precede salon I'invention a 
Taide de moyens de traitement appropries formes par exemple par des proeasseurs de traitement du signal 
ou tout dispositif Equivalent connu est figure sur I'organigramme de la figure 7. 

II consiste suivant T etape 9 a estimer la matrice de congelation Rxx suivant la relation (7). Cette matrice 
4$ est estimee sur un bloc contenant un nombre d'echantiUons K superieur au nombre d'echantillons de la se- 
quence de, synchronisation. 

L'estimation de R^x est utilisee ensuite pour toutes les positions de synchronisation incluses dans le bloc. 
Puis la matrice de correlation inverse R« est cateulee a Tetape 1 0. qui est def ini par la relation (8) est cal- 
cule a r etape 11. Pour chaque echantillon correspondant a une position de synchronisation P dans la suite 
50 des echantillon© P = 1 a P = IC-K, fx* calcule pour i = 1 a M sur M blocs de Ko symboles comma defini a 
la figure 5 suivant la relation 

ss fcai-^- Z X(P+K+iK 0 )d (K + iK 0 ) (46) 

K 0 K =1 
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Des calculs de transfbrmee de Fourier Rapide sont effectues d i'etape survante 12 sur un nombre deter- 
mine M (M = 8 par exemple) de vecteurs suivant la relation (41). 

Aux etapes 13 et 14 un calcul de vecteur R Xfl | decali et un calcul de transformee de Fourier rapide suivant 
la relation (43) sont executes pour ameliorer la resolution par un facteur de 2. 
$ Le critere de synchronisation est calculi a I'etape 15 en application de la relation (44). 

La comparison du critere ou seuil 11 est effectuee a I'itape 16. Si a Tetape 16 la valeurdu critere 




est superieure ou egale au seuil, la position de synchronisation est difinie i I'etape 17 sur la position P de 
T^chantillon correapondant ^estimation du dicalage Ooppler peut ensuite etre aff ini & I'etape 1 8 en calculant 
a nouveau la valeurdu critere de synchronisation pour des dicalages Doppler entourant le decalage 



15 




20 pricident Cette itape peut etre suivie par une itape 19 d'estimation de la position de synchronisation pour 
les autres trajets. 

La recherche de synchronisation priced emment decrite est effectuee pour chacun des echantillons de 
signal numerise. En pratique, pour amiliorer la precision de la position de la synchronisation il est preferable 
d'effectuer un surichantillonnage du signal d'un facteur determine U L itant egal a 4 par exemple. Les trai- 
ls tements de rcrganigramme de la figure 7 et dicrits ci-dessus sont encore utilisables sur le signal ainsi suri- 
chantillonni, mats il est possible d* uttUserune variant© dans le proced&, afin de riduire la puissance de calcul 
nicessaire au traitement Dans ces conditions, aprta avoir calculi pour chaque ichantillon le vecteur r X di pour 
i = 1 f ... M sur M blocs de Ko symboles suivant la relation (46) un critere de synchronisation C t est calculi sur 
chacun des M blocs par la relation 
30 C, = f^R^fx,,, (47) 

Puis un critere somme de i ensemble des est calculi suivant la relation : 



C 1 (P) = XC i (48) 

i=1 

Ce critere correspond a une integration non coherente des Ci (contrairement d ce qui est fait pour le calcul 
par FFT oti on intigre de facon coherente), ce qui se traduit par une remontee du entire sur les echantillons 
-*o de bruit 

Le calcul du entire a lieu par blocs de L echantillons consecutifs. Sur chacun de ces blocs, I'dchantillon 
qui donne la valeur maximum de C^P) est selectionni pour effectuer le calcul par transformee de Fourier ra- 
pide necessaire a I'estimation du d teal age Doppler comme dija dicrit aux etapes 11a 19 de la figure 7. 

45 

Revindications 

1 . Procedi permettant & un modem de se synchroniser sur un transmetteur de signaux nurneriques par voie 
hertzienne en presence de brouilleurs eonsistant a : 
50 . estimar (9)Ja matrice de corrilation R** de signaux regus sur un ensemble de N capteurs de recep- 

tion (1 t ... 1m) par Woes composes d'un nombre d'ichantillons d'une sequence de synchronisation 
transmise par le transmetteur, A 

- a caiculer (10) la matrice de correlation inverse R«, 

- a caiculer (11) des vecteurs d'intercorrelation entre les signaux X(k) recus sur I'ensemble des N 
55 capteurs et un signal de replique d(k) connu. 

- a caiculer un crit&re de synchronisation (7) effectuant !e produitscalaire entre le vecteur d'intercorre- 
lation r w et le vecteur obtenu par produit de ('inverse de la matrice de correlation R^J et du vecteur 
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- & comparer (8) la valeur du critere obtenu a une valeur de seuil ti determinee pour placer la synchro- 
nisation sur Tectiantillon pourlequel la valeur du entire depasse la valeur de seuil ti, caracterise en 
ce qu'il consiste pour estimer le decalage Doppler inherent aux communications entre emetteurs et 

5 recepteurs, a decouper la sequence de synchronisation en un nombre determine de M tranches de 

ko symboles, a effectuer (11 « 1 3) sur ces tranches des calculs de vecteurs d'intercorrelation partielies 
rx« (i s 1. M) de telle sorte que sur chacune des tranches de ko symboles, le decalage Doppler 
fasse toumer la phase d'une valeur inferieure & 90°. d effectuer (12, 14) des calculs de transformee 
de Fourier rapide (FFT) sur les vecteurs ?xdi obtenus, et a calcuier une valeur du critere de synchro- 

10 nisation a partir des vecteurs rw(kAf) A obtenus en resultat des calculs de transformee de Fourier ra- 

pide, et de la matrice de correlation R** 

2. Precede selon la revendication 1 caracterise en ce que Testimatton de la matrice de correlation est 
obtenu par Sa relation 



15 



20 



25 



et le vecteur d'intercorrelation est estirne sulvant la relation 

K 



1 



30 



ou p. est la position du premier echantillon a partir duquel la correlation est effectuee entre le vecteur X 
et la replique d. 

3. Precede selon la revendication 2 caracterise en ce qu'il consiste a choisir la case de calcul de la trans- 
formee rapide qui donne une valeur maximum du critere de synchronisation et & comparer la valeur maxi- 
mum obtenue a une valeur de seuil determines pour decider de la synchronisation. 

35 4. Precede selon Tune quelconque des revendications 1 a 3 caracterise en ce qu'il consiste pour obtenir une 
meilleure resolution surle decalage Doppler et une valeur plus importance du maximum du critere, a cal- 
cuier (13) des vecteurs r x<a decales en frequence eta effectuer un calcul de FFT sur les vecteurs r^ d lo- 
cales de facon a obtenir la valeur du critere avec une meilleure resolution sur les frequences Doppler. 

40 5. Precede selon Tune quelconque des revendications 1 d 4 caracterise en ce qu'il consiste a surechantil- 
lonner les signaux recus sur les N capteurs par rapport au rythme symbol e, pour ameliorer la resolution 
sur la position de synchronisation. 

6. Precede selon Tune quelconque des revendications 1 & 5 caractertsee eh ce qu'il consiste & effectuer le 
45 calcul de la matrice de correlation Rxx sur un nombre d'echantillons K superieur au nombre K d'echan- 

tillons de la sequence de synchronisation af in de diminuer la puissance de calcul. 

7. Precede selon Tune quelconque des revendications 1 & 6 caracterise en ce qu'il consiste pour diminuer 
la puissance de calcul et apres avoir surechantillonne d'un facteur L par rapport au rythme symbols, a 

SO effectuer une etape de pretratternent consistant d calcuier pour cnaque echantillon p (p = 1 a L) composant 

un symbole les valeurs C, du critere de synchronisation sur chacune des M tranches de ko symboles et 
a calcuier la valeur du critere en effect uant la somme des cri teres C, obtenus pour retenir dans la suite 
du traitement Cechantillon donnant la valeur maximale du critere C - ZC». 

35 8. Dispositif pour la mise en oeuvre du procede selon I'une quelconque des revendications 1 a 7 du type 
comportant un reseau de N capteurs (1 ,...1 N ) de signaux de reception couple a un reseau de f iltres adap- 
tes (2) comportant des circuits de ponderations amplitude phase comma n des par un processeur de 
commande (3). le processeur de commande (3) etant programme pour 
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- estimer la matrice de correlation inverse R£ des signaux repu et calculer les vecteurs cTinter- 
corr&ation fxa entre les signaux reQus et le signal de r$plique d^ caract6ris6 en ce qu'il comprend 
des moyens de traitemem (5...19) pour calculer un critere de synchronisation (7) en effectuant un 
produit scalaire entre le vecteur (Tintercoirelation et le vecteur obtenu par le produit de llnverse 
de matrice de correlation R« et du vecteur d'intercorrelation 

- comparer (8) la valeur du entire obtenu a la valeur de seuil i\ pour placer la synchronisation sur 
rfechantillon pour 1 equal la valeur du critere de passe la valeur du seuil r|, 

- et estimer le decalage Doppler inherent aux communications entre 6metteurs et recepteurs. en de- 
coupant la sequence de synchronisation en un nombre determine de M tranches de ko symboles, 
en effectuant (11, 13) sur ces tranches des calculs de vecteurs d^intercorr^ation partial es 0 = 

1 M) de telle sorte que sur chacune des tranches de Ko symboles le ddcalage Doppler fasse tour- 

ner d la phase <fune valeur irtterieure & 90°, en effectuant (12, 14) des calculs de transforms de 
Fourier rapide (FFT) sur les vecteurs obtenus, et en calculant une valeur du critere de synchro- 
nisation £ partir des vecteurs rxd(r%£f^ obtenus en r£sultat des calculs de transform6e de Fourier ra- 
pide et de la matrice de correlation R*x- 
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